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ABSTRACT 

This  report  is  concerned  with  the  incompressible, 
Irrotational  flow  of  fluid  in  circular  channels  or  circular 
troughs.   Fluid  circulates  through  these  channels  and  is 
contained  in  them  by  the  action  of  a  gravitational  body 
force.   Under  certain  velocity  conditions,  permanent  waves 
of  large  magnitude  and  constant  angular  velocity  may  be 
solutions  to  the  differential  equations  of  motion.   These 
waves  are  closely  related  to  the  cnoidal  and  solitary  waves 
of  the  infinite  straight  channel.   Expressions  for  their 
shape,  the  critical  velocity  condition  and  their  relation 
to  waves  in  straight  channels  are  derived  and  discussed. 

The  irrotational  case  is  treated  thoroughly  and  the 
equations  for  the  rotational  case  are  derived  and  one 
solution  is  developed  in  detail. 
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Chapter  I 
INTRODUCTION  -  CONCLUSIONS 

Since  l844  when  Scott  Russel  made  his  first  observation  of 
a  solitary  wave  on  a  stream,  much  theoretical  work  has  been  done 
on  the  problem  of  solitary  waves  in  a  liquid  of  constant  density. 
Almost  all  of  this  work  has  been  in  two  dimensional  flow,  or  in 
Infinite  straight  channels.   There  has  also  been  some  work  done 
on  the  associated  cnoidal  waves  for  these  surfaces  and  also  on 
the  surfaces  of  infinite  cylinders  and  toroidal  channels  with 
radial  gravitation.   In  this  report  we  shall  work  with  a  different 
geometry  (the  circular  channel  or  bowl)  which  is  described  later  in 
this  chapter. 

In  general,  the  solutions  will  have  several  things  in 
common  with  previous  results.   The  first  of  these  is  the  charac- 
teristic critical  speed.   For  certain  velocity  profiles  the  equa- 
tions will  exhibit  a  bifurcation  i.e.,  there  will  be  two  families 
of  solutions.   The  first  family,  which  will  be  called  the  steady- 
state  flow,  will  exhibit  no  wave  in  the  direction  of  the  flow. 
This  family  of  solutions  is  also  common  to  flows  at  non-critical 
speeds  and  is  an  exact  solution  of  the  differential  equations  of 
fluid  flow.   The  second  family  of  solutions  will  be  waves  of  a 
permanent  type,  that  is,  the  basic  shape  of  the  wave  is  fixed 
although  it  may  move  over  the  surface  of  the  water.   The  second 
thing  in  common  with  the  previous  work   in  this  subject  will  be 
the  shape  of  the  wave;  it  will  in  general  be  the  characteristic 


cnoldal  profile  (Jacob!  Elliptic  Function)  In  the  direction  of 
the  flow.   A  third  thing  In  common  with  previous  work  is  that 
the  specific  wave  shape  depends  on  the  solution  of  an  elliptic 
Neumann  problem  in  the  equilibrium  or  steady  state  cross  section 
of  the  fluid. 

The  steadily  progressing  waves  discussed  here  will  be  waves 
whose  shape  remains  fixed  as  viewed  from  a  coordinate  system  which 
is  rotating  with  respect  to  the  fixed  coordinate  system  at  a  con- 
stant angular  velocity.   A  special  case  of  this  will  be  stationary 
waves,  i.e.,  where  this  angular  velocity  is  zero.   Let  me  note 
here  that  we  will  depart  from  the  main  body  of  previous  works  by 
describing  all  these  waves  from  a  fixed  coordinate  system.   We 
will  find  this  to  our  advantage  within  the  confines  of  the  geometry 
of  this  problem. 

We  will  be  interested  in  waves  in  circular  channels,  see 
Diagram  2-1.   Consider  a  fixed  curve  (s-G(p)  =  0)  in  a  cylindri- 
cal coordinate  system  where  (p)  is  the  radial  coordinate,  (a)  is 
the  angular  coordinate  and  (s)  is  the  axial  coordinate.   This  fixed 
curve  is  rotated  about  the  axis  of  the  coordinate  system  and  thus 
it  generates  a  trough  or  circular  channel.   Fluid  will  be  contained 
in  this  trough  by  the  action  of  an  axial  constant  gravitational 
force. 

It  is  not  a  priori  clear  that  cnoidal/solitary  waves  exist 
in  circular  channels  since  their  existence  in  straight  infinite 
channels  is  strongly  coupled  with  the  fact  that  waves  of  very  large 
wave  length  exist  in  such  channels.   The  infinite  length  of  the 


straight  channel.  In  which  It  has  been  proved  that  cnoldal  waves 
can  exist,  stands  In  contrast  to  the  finite  length  of  the  liquid 
filled  circular  channel  considered  here.   The  existence  of  these 
waves  In  finite  channels  Is  suggested  by  observations.   In  par- 
ticular, the  solitary  waves  observed  for  finite  times  under 
natural  circumstances  and  also  under  laboratory  conditions  in 
finite  channels  strongly  suggest  that  these  permanent  waves  do 
exist.   It  is  for  these  reasons  that  we  apply  slight  variations  of 
the  techniques  for  infinite  straight  channels  to  the  problem  of 
circular  troughs.   What  we  do,  basically,  is  associate  a  stretching 
parameter  to  some  of  the  dependent  and  independent  variables  so 
as  to  emphasize  the  angular  direction  and  angular  velocity.   We 
then  expand  each  of  the  new  dependent  variables  in  a  power  series 
in  this  parameter.   The  basic  assumption  is  that  these  power 
series  converge  to  solutions.   We  substitute  these  series  into  the 
equations  of  fluid  motion  and  then  extract  the  equations  coeffi- 
cient to  different  powers  of  the  stretching  parameter.   That  is, 
the  equation  L(x,  e,(|){x,  e )  )  =  0,  where  (L)  is  a  differential  opera- 
tor (or  boundary  condition),  and  (x)  represents  time  and  space 
variables,  and  (e)  is  the  stretching  parameter,  and  (<j))  is  the 
desired  solution  will  yield  the  following  expression  under  the 

assumption  ^{x,e)   =  y        ^. (x)e  . 

1=0 

oo  ^ 

L  =  y~  L.  (x,(t)  ,  .  .  .  ,(|)  )e   =0 
1=0 

In  this  expression  the  (L. ).are  differential  operators  (or 


boundary  conditions)  Independent  of  (e).  These  resulting  equa- 
tions are  solved  recursively  starting  with  those  coefficient  to 
the  lowest  powers  of  the  stretching  parameter.   This  assures  us 
of  consistent  results  from  all  the  equations.   The  results  in  all 
cases  considered  show  that  the  leading  terms  in  the  expansion  of 
the  wave  profile  have  a  marked  cnoidal  shqpe  in  the  angular  direc- 
tion.  The  shape  of  the  wave  in  the  radial  direction  Is  highly 
dependent  on  the  assumed  steady  state  fluid  velocity  and  the 
velocity  of  the  progressing  wave. 


Chapter  II 

EQUATIONS  OP  MOTION: 
ROTATIONAL  AND  IRROTATIONAL 


Let  the  spacial  coordinates  be  {p,a,s)  and  let  them  be  Indexed 
(1,2,3).   The  velocity  vector  (V)  in  this  cylindrical  coordinate 
system  is 

( 2-1 )  V  =  T^V^  +  TgVg  +T^V^  . 

/o  o\  T7-    dp     ,j  da  ds 

(2-2)  V^  =  dt  '    ^2  =  Pdt  '    ^3  =  dt  • 

The  velocity  (V)  will  be  referred  to  in  the  body  of  this 
report  as  the  fluid  velocity  and  as  the  particle  velocity. 
In  this  coordinate  system,  we  have 

(2-3)  ^  =  ^1^-"^^}^^^^^  ' 

t2-M  V.V=i^(pV,)  +i^(V2)  ^-^(V^)  , 

(2-5)   vxv=T^civ  -v^I+T^Cv^  _v  }+i  [l^(pV2)-iv  )  . 

f^^a    s        sp    ^r-r       "^a 

The  equations  of  the  free  surface  and  the  fixed  surface  will 
be  given  respectively  by  , 

{2-6)  ■  s  =  Ti(p,a,t)  , 

(2-7)  s  =  G(p)  . 


The  equations  of  continuity  and  momentum  are  given  respec- 
tively as 

(2-8)  V.  V  =  0  , 

(2-9)  F  -  ^   gFid  P  =  ^  . 

In  these  equations  (P)  is  the  pressure,  (t)  is  time,  (6)  is 
the  constant  density  and  (F)  is  the  body  force.   To  be  specific, 
we  will  solve  the  problem  with 

(2-10)  F  =  -T^g  , 

where  (g)  is  a  positive  constant.   We  can  think  of  the  fluid  held 
in  a  circular  channel  (eq.  (2-7))  by  the  action  of  a  gravity  force 
(eq.  (2-10))  directed  along  the  axis  of  symmetry  (p  =  0).   See 
Diagram  2-1. 

In  (2-9)  we  complete  the  differentiation  on  the  right  using 
relations  (2-2);  this  together  with  (2-8)  gives  the  set  of  lOur 
simultaneous  differential  equations  which  we  will  solve. 

(2-11)  -  |-(pV,  )    +l^v^+^V,  =  0 

^  '  p   dp   '^    1"^         p   da      2        ds      3 

(2-12)  V       +V,V      +—  VV       +VV       -i(V)=-ip 

W  P  Ob  S 

f2-l'5)  V       +VV       +—  VV       +VV       +  —(V   V    )    =    -   ^  P 

\^   i:>)  ^2.         12        p   ^2^2    ^  ^32        P^    1   2^  ^  ^a 

up  ot  s 


f2-l41  V       +V  V       +—  V  V       +V  V  =    -  —  ?     -  a 


We  will  solve  these  partial  differential  equations  subject 
to  zero  pressure  on  the  free  surface  and  the  natural  kinetic 
conditions  (2-15,  l6,  17).   On  s  =  Ti(p,a,t), 

(2-15)  P  =  0 

and 

(2-l6)  Y      =  T\   Y^  +-  r]   V„  +Ti^  , 

3    'p  1   p  'a  2   't  ' 

and  on  s  =  G(p), 

(2-17)  V,  =  ^^  V,  . 

3    dp    1 

To  reiterate,  the  problem  Is  to  solve  the  set  of  relations 
(2-11)  through  (2-17). 

We  now  introduce  a  stretching  technique  to  facilitate  finding 
a  solution  to  the  problem.   Coordinates  (p,a,s)  go  into  (p,9,s), 
and  (P,V^,V2,V^,t)  go  into  (TT.U,V,W,t)  respectively. 

(2-18)  9  =  /^a 

T  =  /eg  t 

U  =  V^/i/g 
V  =  Y^T/g 

W  =  Vy'eTg 

TT  =  P/5g 

With  this  change  of  variables,  the  problem  is  now  to  solve 
the  following  relations  simultaneously,  and  then  to  use  (2-l8)  to 
identify  the  variables. 


(2-19)     e^(V^)  +^(pU)p  +  W3  =  0 

(2-20)     e[U^  +  Tf^  +  i  VUg  -  ^  V^]  +  UUp  +  WUg  =  0 
(2-21)      e[V^  +  i  7]-0  +  ^  Wg]  +  iuv  +  UVp  +  WVg  =  0 
(2-22)     e[W^  +  TTg  +  ^  VWg  +  1]  +  UWp  +  WWg  =  0 
On  the  free  surface  {s  -r]{p,9 ,t)   =  0]    we  have 

(2-23)  TT  =  0  , 

and 

(2-24)  [W   -   UTip]    =    e[Ti^   +  VTig/p]    . 

On  the   fixed   surface    (s-G(p)    =  O)    we  have 

(2-25)  W  -UG'(p)    =   0    . 

[Note  for  the  body  of  this  report  we  will  use  the  following 
alternate  conventions:   ^4-^=  ('),  ^4^=  (')•] 

We  also  note  that  the  vorticity  in  these  variables  is 

(2-26)   V  X  V  =  /i/I   [/?  T^(i  Wg-Vg)  +T2(U3-Wp) 

+  /rT^(^^(pV)  -iu,)]  . 

Of  special  Importance  is  the  case  of  irrotational  motion. 
The  relation  (V  x  V  =  0)  implies  that  (U,V,W)  can  be  written  as 
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(2-27)  U  =  Jp 

Where  []^  =  ^(p,9,s,T)]  is  the  velocity  potential.   Equations  (2-19) 
through  (2-25)  now  become 

(2-28)        ^l^^.i^(p|j)l.l^^   =   0 

r 


2. 
-pr    '    '  '  P    '   "p^  ^e  ^Sp   "   p3  ^9    '    '  ^p  ^pp    '   ^s  ^sp 


(2-29)         e[i       +TT     +-iTLlo      -^Inl+$     I       +!.$       =0 


(2-31)        £[?       +  TT     +  ^  ?.  "^^     +1]   +??       +?'^       =0 

r 

Ons   =  Ti(p,0,T), 

(2-32)  TT  =  0    . 

On   s  =  T\{p,9,T) , 

(2-33)  tls-5p^p^   -■=   ^[^+^^0^0^ 

^   ^  P 

On   s   =  G(p), 

(2-34)  Is-$pG'    =  0    • 


We   Integrate   equations    (2-29,    30,    31)    to   get   the   energy 
integral    (2-35). 


(2-35)     ^[5,  -  TT  -  s  .  -1^  1,2]  .  1  f J^2  ^  1^2^  ^   E(,j  . 

This  simplifies  the  problem  of  solving  this  set  of  equations  for 

:•:,  T.   We  now  must  solve  the  set  (2-28,  33,  34,  36)  and  then 
the  pressure  is  determined  from  (2-55). 


(2-28)  4v^..  +  i^  D  ^ 

On  s  =  n(p,0,T), 


75ee-^^P^$-f3s  =  ° 


On  s  =  ti(p,0,t). 


(2-36)    [E  -  ^(f^)2  _  i(j  )23  ^  ,f^  ^  J  ^  -^(le)'] 

^  2p 

On  s  =  G, 


(2-34)  L-g'I  =  0  . 
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Chapter  III 
IRROTATIONAL  STATIONARY  CNOIDAL  WAVES 

In  this  chapter  we  will  look  for  solutions  of  the  partial 
differential  equations  which  are  time  independent.   Let 


(5-1) 


1,  -  ^r  =  \  =  ^    ' 


then  (2-28,  33,  3^,  35,  36)  become  (3-2,  3,  ^,  5,  6) 


(3-2) 


On  s  =  Ti(p,0), 


(3-3) 


7^   P^'^   ^ 


^  =  0 


^  ^   p 


and 


(3-4) 

and  on  s  =  G(p) , 

(3-5) 


ri  =|[E-^($3)^  -H)']  -^^le)'  ' 


Is  -  Ip^'  =  ° 


Equations  (3-2,  3,  4,  5)  are  to  be  solved;  then  the  pressure 
is  derived  from  (3-6). 


(3-6) 


^  =  i[E-^(fp)"-i(S,)2l  -^($9)"-s 
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Before  proceeding  with  the  solution,  let  us  Identify  some 
constants  and  symbols  with  respect  to  Diagram  5-1. 

Let  (a,  a+b)  be  defined  as  the  abscissae  associated  with 
the  intersection  of  the  curves  (s  +  7  /2p   =  0)  and  (s-G(p))  =  0). 
Let  (  ^)  be  the  surface  contained  between  those  two  curves.   (S^) 
will  denote  the  boundary  of  (  ^ ) .   It  will  be  shown  later  that  (  <£) 
will  be  the  cross  section  of  the  trough  which  will  be  occupied  by 
the  steady  state  fluid  rotating  about  the  axis  (p  =  0).   It  will 
be  necessary  later  to  use  the  divergence  theorem  on  this  region; 
the  specific  form  is  derived  here. 

On  (s  +7  /2p  =0),  the  unit  outward  normal  (N)  and  the 
differential  length  (djrf)  along  the  boundary  [h  Si)   are 


(3-7)  N  =  (^ 


P^ 


Ixf  =  ]l  +  ^  (-dp)  . 


d;J  =  J 

P' 

On  (s-G  =  0),  (N)  and  (d^f)  are  given  by  (3-8) 


(3-8)  N  =  (g',  -l)j  ]i+  (G')^ 

di/=  il+  {G'f    (dp) 
The  divergence  theorem  will  now  read  as  (3-9)  or  (3-10) 
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Oio)  ff 


L^^p^p)  +^^P^s^ 


dA 


a+b 


p(Xc 


2 
7   X. 


a+b 


dp 


/   P(>^s-^'Xp) 


~7 


dp 


s=G 


2p 


We  will  now  proceed  with  the  formal  solution  to  equations 
(3-2,  3,  4,  5).   Let 


(3-11) 


(3-12) 


(3-13) 


f(p,e,s;e) 

^(p,Q;e) 


OD 


^$  (p,0,s)e^  , 
i^  ^ 


00 

ZZ  Ti.(p,0)e^  , 


OD 

i^ 

00 


E  = 


i=0 


where  the  (E. )  are  constants. 


To  begin,  we  will  examine  the  equations  whose  coefficient 


0 


is  (e  ).   Equations  (3-2,  3,  ^,    5)  imply 


(3-14) 

on  s  =  G, 

(3-15) 

and  on  s  =  r]„. 


f^^P^)  ^^fP^)]^  =  0 


5q  -^   G'(p)  =  0  , 
s    p 


Note:   We  recall  that  the  following  alternate  conventions  shall 

^(  )  =  (•),  ^^ 


be  used:   ^(  )  =  ( *  ) ,  :^(  )  =  (') . 
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(3-16)  lo   -%%    =  0 

s         p     p 


and 

2 

P 


(3-17)  Eo=^[lo2+V 


The  first  three  of  these  equations  imply  that 

(3-18)  %{p.Q,s)   =   ^(9)  . 

That  is,  {^)  is  a  solution  of  a  strongly  elliptic  equation 
in  region  ( /?  )  with  zero  normal  derivative  on  the  boundary  (S^). 
The  maximum  modulus  theorem  applies  and  we  have  (3-l8).   For  con- 
venience we  will  rename  (^)  in  (5-19).   Equation  (3-17)  now  gives 
us  (Eq  =  0). 

(3-19)  lo  =  ^0^^^  • 

Let  us  now  examine  the  (e  ) -equations.   Equations  (3-2,  3,  ^, 
5,  19)  imply 

on  s  =  G, 


(3-21) 

on   3   = 

^0' 

(5-22) 

and 

II  -$1  G'  =  0  > 
s    p 


5l  -II  ^0  =  Fq^o  /2P^  ' 
S     p   p        a 
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(5-23)  ^0  =  ^1  -  -^(^o)^  • 

2p 

We  differentiate  (3-23)  and  substitute  It  In  (3-22)  and 
get  (3-24).   On  s  =  t^q, 

3 

P  "P        2p 


(>24)         Ii  -Ii  no_  =  _ro_^^o 


We  use  the  divergence  theorem  (3-10),  identifying  (x)  with 
i^-,)   and  conclude  that  (F„  =  0).   Therefore  (3-22)  can  be  restated 
as 

(3-25)  "^0  ^  %  "  ^^/^P^  ' 

where 

(5-26)  t)  =  ^0  =  ^^  • 

To  describe  all  the  solutions  for  a  specific  channel  (s  =  G) 
It  is  necessary  to  keep  the  constant  (E-,  )  arbitrary   throughout 
all  the  succeeding  equations.   This  is  cumbersome,  and  since  we 
are  interested  in  all  the  solutions  for  all  possible  shapes 
(s-G  =  0),  it  is  somewhat  redundant  to  keep  (E-,  )  arbitrary.   Let 
us  set  (E-|  =0).   In  doing  this,  we  must  remember  that  the  group 
of  solutions  for  a  specific  channel  (s-G  =  0)  will  be  replaced  by 
the  totality  of  solutions  for  the  group  of  channels  (s-G  -A  =  0), 
where  (A)  is  an  arbitrary  constant. 

We  now  have 

(3-27)  tIq  =  -7^/2p^  . 
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The  equations  (3-20,  21,  22)  for  ( J-^ )  are  identical  to 
(3-14,  15,  l6)  for  (^)«   For  the  same  reasons  as  before  we  get 

(3-28)  f^(p,9,s)  =  f^iO)    =   F^(e)  . 

Let  us  now  examine  the  (e  ) -equations.   Equations  (3-2,  3, 
4,  5)  Imply 

On  s  =  G, 

(3-30)  $P  -  Jp  G'  =  0  . 

s    p 


On   s  =  T^Q  = 

-7^ 
2P^ 

(3-31) 

and 

(3-J2) 

$2  -^2  ^0   =  (^^1«^/P^  ' 
s     p   p       9 


Til  =  E^  -  ^^l/P^ 


We  combine  the  last  two  equations  to  get  on  s  =  r]    , 

(3-33)  I2  -^2  "^0  =  -y^i/p^   ' 

s    p  p 

Using  the  divergence  theorem  (3-10)  on  equations  (3-29»  30, 
33)  we  get 

a+b   >^ 
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Equation  (3-34)  represents  a  bifurcation  typical  to  this 
problem  and  will  reappear  in  subsequent  chapters  in  a  similar 
form.   If  (F-,  =0),  we  get  irrotational  flow  with  no  wave,  i.e., 
the  trivial  steady  state  solution.   That  is 


n  =  TIq  "^  ^^1  ^  °^^  ^ 


-7 


+  e 


77- 


E. 


+  0(6^) 


=  T){p)  +  0(e'-)  , 


where  7-,  =  F,  =  const.   If  similar  assumptions  are  made  throughout 
this  expansion  we  get 


Tl^  =  n^Cp)  • 


The  final  result  is  (tj  =  T](e,p)  =  A  +  B/p  )  where  tA,B)  are  con- 
stants with  regard  to  the  space  and  time  variables.   This  is  a 
steady  state  solution  and  satisfies  the  exact  equations  of  motion. 
In  fact  it  can  also  be  represented  as  r]  =  t^  (p).   For  the  simplest 
form  of  the  solution  for  which  there  may  be  a  wave  we  must  look  at 
the  alternate  solution  to  (3-3^).   Assume  (F,  ^  0),  then  (7  )  is 
determined  from 


(3-35) 


dA 


=  7 


/e 


a+b 
a    ^ 


Integrating  with  respect  to  (s)  we  get 


17 


a+b  2 


(5-J6)  fff'f       (-^->P- 

■4     P       J  2p5       P 

Equations   (3-35,    36)    imply 

a+b  /     a+b 

(3-37)  y^  =  -^J       (G/p)dp  //       (l/p^)dp    , 

a  /a 

a+b 


"-38)     ^'=<r)i^^/  t-o/p)^p- 


This  value  of  (7),  which  we  shall  call  the  critical  value, 
can  be  related  to  the  critical  speed  for  solitary  waves  in 
infinite  straight  troughs.   This  can  be  seen  by  letting  (a)  get 
large  with  (b)  relatively  fixed.   Let  (Ht)  be  the  average  height 
of  the  fluid.   In  (3-25)  we  get 


(3-39) 

JJf^^-Ml                 as      ,a,t. 

n 

and 

a+b 

(3-40) 

J      (l/p^)dp  ^b/a^     as      {a)t    • 
a 

Therefore 

7^  —  a^(Ht)                       as      (a)t    . 

The  tangential  velocity  squared  is  given  by 


K' 


(3-41)  v^-f-^j  =h-*h   ^^    (^)^ 
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This  is  seen  to  mean 

(3-42)  V^  ^  (Ht)   as   (a)T 

In  dimensional  form,  this  Is 


(3-^3)  (l^-  V)  -/(gHHt)   as   {a)t 

This  is  the  same  as  the  critical  speed  for  solitary  and  cnoidal 
waves  in  infinite  straight  channels;  hence  the  term  critical  value 
of  (7). 

Restricting  the  value  of  (7)  to  the  critical  value  given  by 
(3-38)  we  can  now  proceed  with  the  solution  of  (3-29,  30,  32,  33). 
Using  the  maximum  modulus  theorem,  we  see  that  (^p)  can  be  written 
as 

(3-44)  $2  =  F^(9)f(p,s)  +  F2(0)  , 

where 


(5-^5)  [^p^.|_p^' 


^  =  -1/p  , 


and  on  s  =  G, 


(3-46)  Y  -  Y  g'  =  0 

'  s    p 


and  on  s  =  x]^, 

(3-47)  ^3  -  ^pTiQ   =  -y^/p"^ 


with  (Fp(9))  arbitrary  up  to  this  point. 
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The  solution  of  (3-45,  46,  4?)  is  not  obvious,  even  for 
specific  troughs  (G  known  explicitly).   If  we  solve  (3-45,  46,  47) 
approximately  or  exactly,  we  will  be  able  to  solve  for  (F, (9)) 
approximately  or  exactly.   We  can  determine  the  nature  of  (F,(0)) 
without  knowing  (^)  explicitly;  for  this,  it  is  necessary  to  go  to 
the  next  higher  order  equations  (coefficients  of  (e-^)).   These 
equations  are 

(J)  ••••         •• 

(5-^8)   (^  p  ^  +  ^  p  ^)5,  =  .  -^  =  .  y  -  !|  . 
On  s  =  G, 


(3-49) 

On  s  =  tIq, 

(3-50) 


$   -  $3  g'  =  0  . 
^s    ^p 


2/^  2. 


TI2  =  E^  -  [(F^)V2p^]  -rl^     /?      , 


and 


051 )  %       -    TIqI, 
s     p  ^p 


7F^F^ 


1. 


ss 


7^Y 


£S_ 


L  P 


^r  ^  3 
p    p  . 


p   ....       Oj, 


^  +  V2t%s^0   -  ^ss^ 


p  p  ^^     p 


To  get  equation  (3-51)  we  used  (3-3)  together  with  the 
following  relations. 


(3-52) 


i 


-\ 


+  e^l$S 


ss 


+  O(e^) 


'0 
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053) 

On  s  =  tIq, 


fc 


=  % 


+ 


£^1$. 


'0 


-ps 


+  0(£^) 


'0 


^7      2 

^   -  -^  =  ^  • 
S      3       4 

P       P 


We  now  apply  the  divergence  theorem  to  (3-48,  49,  51)  and 
get  (3-55). 


(5-55)   F;W/IdA-r 


a+b 
2  r    Y 


dp   +F^     ^-7 


S=T) 


0 


2/JJ   p 


a+b 
dA   .2  f    d£ 


.^^^  (2p-2p^%+7^-yp3-P^^ss) 


^^AfJ     IF 


dp 


s=ri 


0 


a+b   3„,    2,„ 
n  p  T   -7  Y 

+  ^1^2/,/     — ^V — ^ 


dpf   . 


S=T] 


0 


2 

Since  (7  )  is  the  critical  value  (3-35)  we  see  that  the 

coefficient  of  (Fp)  in  (3-55)  is  zero.   Let  us  now  examine  the 
coefficient  of  (F, ). 


(3-56) 


-  dA  = 
P 


Y(V-  pW)dA 


^ 


/? 
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(3-57) 


a+b 


S=T1. 


a+b 
a 


s=ii 


(-  \)dp 
P 


a+b 


=  /    (P^)(v 


dp 


a=^, 


P^^N^ 


-^ 


=  //  (V.  {p^W)dA 


Combining  (3-56,  57)..  we  t^et  the  following  result 


(3-53)      rri« . 


Let 


a+b 
2    f  T 


a 


7 


dp 


S=T1, 


=  //  (V  .  (p^W)  -  y(v.  pW)) 


dA 


^ 


=  /7  p(W-  W)dA 

=  ^p[(^p)^  +  (^3)^]dA  >  0  . 


(3-59) 
(3-60) 


Y{9)    =  F^{9) 


m, 


0 


=  ^p[(^p)^  +  (^g)^]dA  >  0 


a 
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(3-61) 


(3-62) 


m^  =  7 


r^'^C   2p-2p^^  +7^¥   -p^Y 
'    »   ^  ^      p   _  ps  ^   ss 


a+b 


nig  =  Eg 


p^ 


7^^ 


ss    p2 


£^ 


s=t1q=-7  /2p 


dp 


dp 


s-r\ 


0 


Using  these  definitions  and  the  previous  results  we  find 
that  equation  (3-55)  can  be  rewritten  as  (3-63)* 


(3-63) 


hIqY  =  m  YY  +  nigY 


Equation  (3-63)  occurs  in  previous  works  on  solitary  waves. 
It  will  also  appear  again  in  subsequent  chapters.   For  these 
reasons  we  will  discuss  in  detail  the  applicable  solutions  of  this 
equation  for  our  problem.   Let  us  first  reiterate  some  details. 
For  a  specific  shaped  trough  we  must  determine  (a,b,7  )  according 
to  Diagram  3-1  and  equation  (3-35)  o^  (3-38).   We  then  solve  the 
elliptic  Neumann  problem  (3-^5,  46,  47)  for  (Y).   Coefficients 
(mp^,  m^  ,  mp)  are  determined  from  (3-6o,  6l,  62)  once  (^)  is  known. 
We  then  can  solve  (3-63)  for  (Y(0)).   The  solutions  of  (3-63)  fall 
into  three  categories  depending  on  the  value  of  the  coefficient  of 
the  nonlinear  term  (m-,). 


Case  1:   m. 


0 


Equation    (3-63)    now  reads   as    (3-64). 


(3-64) 


itIqY  =  mgY 
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ly  smooth  solutions  occur  If  (nu  "*  0);  then  we  have 
sinusoidal  waves. 


r,  +  r 


r^  -  r. 


(-  '  •     Y  =  (  ^g  ^)  +  (  ^g  ^)  cos 


m, 


m. 


0 


(0  -9*) 


There  is  a  restriction  on  (nu/mQ)  so  that  (Y)  will  be  a  continuous 
solution,  the  restriction  implies 


{3-66) 


-(J)  m^  =  m^e  ;     (J)  =  1,2,3,..., 


(5-67) 


-(J)' 


p[(^3)^  +  (^p)^]dA 


a+b 


=  eE 


2 

7    "H 


P^ 


£S 


ss 


2 
P  J 


dp  <  0 


S=T] 


0 


This  restriction  now  applied  to  (5-65)  gives  us  (3-68) 


(3-68) 


Y  =  (^^^^)  +  C^'/^)    cos  (J)a 


A  final  restriction  based  on  the  volume  of  fluid  being  the  same 
with  this  wave  as  without  it,  gives  us  (3-69). 


(3-69)         E2(2ab  +  b^)  =  7{r^  +  r2)  In  (l+b/a) 


Combining  these  results  gives  us  this  final  expression  for  the 
leading  terms  of  the  expansion  of  (i^). 
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(3-70)  x]   r^  t]q  +  er]^   =   - 


7 


2p' 


/ 


a+br 


pY 


2 


ss 


p  J 


dp 


s=r| 


0 


2ab  +b 


2p'^  In  {1+  b/a) 


P cos  [(j)(a-a*)] 

2p 


In  conclusion.  If  (m-,  =  0)  we  need  only  specify  (j),  the 
number  of  cycles,  and  [e(rp-  r, ) ] ,  the  total  height  of  the 
disturbance.   This  is  enough  to  determine  ("npi+eT^-,)  unique  up  to 
an  arbitrary  rotation  (a*),  and  is  given  by  (3-70). 

Case  2  and  Case  3  are  almost  Identical.   Case  2  is  charac- 
terized by  (m-,  <  0)  and  Case  3  is  characterized  by  (m,  >•  0).   To 
solve  for  Case  3  we  merely  transform  (+Y)  directly  into  (-Y),  and 
then  equation  (3-63)  falls  into  Case  2.   Therefore^  we  need  only 
discuss  the  solutions  of  Case  2. 


Case  2:   m,  <  0 

To  get  the  simplest  form  of  the  solution,  we  let 


(3-71) 


Z  =  Y  +  m^/m. 


This  transforms  (3-63)  to 


(3-72) 

Let 

(3-73) 


m^Z  =  m-,ZZ  . 


A  =  -3mQ/m-j^  >•  0  , 
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(3-7M 


c  =  -m^/yEgm^  *■  0 


Then  (3-52)  reduces  to 


(3-75) 


^1  =^2 


1  - 


2^ 
7  c 

— ? 

P  J 


7Z 

P 


We  will  now  show  that  Z{9)   exhibits  the  characteristic 
cnoidal  shape,  and  we  will  determine  which  parameters  are  arbitrary. 
We  begin  to  solve  for  (Z)  by  integrating  (3-72)  or  (3-76),  multi- 
plying by  (Z)  and  then  integrating  again. 


(3-76) 


?v  Z  =  -3ZZ 


(3-77) 


A(Z)^  =  -Z^  +  A^Z  +  Aq  =  P(Z) 


Let  us  look  at  the  roots  of  [P(Z)).   (A  )  is  the  product  of 
the  roots.   Looking  at  (P(0)),  we  see  that  {(sgn  A  )  =  ( sgn  A)). 
This  gives  us  our  first  relation 


(3-78) 


1,2,3 
1  r  (Roots  of  P)  >  0  . 


The  absence  of  the  quadratic  term  in  (p(Z))  Implies  the 
second  relation 


(3-79) 


1,2,3 

>     (Roots  of  P)  =  0  . 


These  two  relations  imply  that  two  roots  are  negative,  one  is 
positive.   Let 


(3-80) 


Ql  >■  qg  *"  0  • 
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The  roots  of  [P(Z)}  can  now  be  written  as  follows: 
(3-81)  +q^,    -qg,    ^2"  ^1  • 

Equation  (3-77)  can  be  written  as  (3-82). 

(3-82)        -h{zf   =   ?{Z)    =  (q3^-Z)(q2  +  Z)(q^-  qg+Z) 

2  2      ^ 

=  q-Lq2(qi  -  ^2^  ^  ^^i  '  ^^1^2  '^^2   ^^  "  ^ 

Without  loss  of  generality,  let  us  order  the  negative  roots.   Let 

q-L  -  ^2  ^  ^2 


or 


q^  ^  2q2 


With  the  roots  thus  specified,  {P(Z)}  appears  as  shown  in  Diagram 
3-2. 

We  make  the  substitution  (3-83)  in  (3-82)  to  get  the 
differential  equation  into  the  more  familiar  form  (3-84). 

2  2 

(3-83)  Z  =  q,  cos  O  -  qg  sin  O 

2 

=  -q2+  (q]_  +  q2)  ^°^    ^ 
o  =  o(0) 

f^   R)l^         ^dn^2     4A      ,     ^1  "^  ^2   .  2  ^ 
(>S^)         (dF)   2q^  -  qg  =  1  -  2q^  -  q^  ^^^^  " 

We  simplify  the  notation  with  the  following  substitution  of 
constants. 
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k  = 


2^1  -  ^2 


?  = 


2^1  -  ^2 


Let  us  note  that  (k  )  and  (i)   are  related  as  In  (3-86) 


(v_   - 


l^^+l  =  3A 


Equation  (3-84)  now  reads  as 


(3-87) 


nTK  dn  _  / 
J^  dF   7 


~5 2 — 

1  -  k   sin  O 


This  equation  is  integrated  now  so  as  to  define  the  function  {O(0)] 


(3-88) 


This  defines  O. 


o 


dO 


0  Ji --■''-  ■' 


k   sin  O 


(3-89) 


O  =  am 


(/ 


(9  -0  ),k 


It  is  of  Interest  to  list  the  bounds  of  {k^,^,^k^)  for 


future  reference. 


Table  3-1 


qg  =   0 

0   <   qg  ^   qi/2 

k^   ^   1/2 

1/2   <   k^  <   1 

f      =    2 

2   >   C   >   3/2 

k^e^l 

1   <   k^   <    3/2 
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We  now  let 


(3-90) 


K  =  K(k) 


E  =  E(k) 

be  the  complete  Elliptic  Integrals  of  (k)  of  the  first  and  second 
kinds.   We  also  let 


(>91) 


H  =  eq, 


be  the  dimensioned  height  of  the  wave.   Then 


(3-92) 


a 


He 

TK 


O 


dO 


0  J  1  -  k^  sin^  O 


or 


(3-93) 


n  =  am 


(»/IL  ■")  =  --(«/! 'i^)- 


The  standard  notation  for  the  trigonometric  functions  of  (O) 
is  given  below. 


(3-9^) 


cn(A,B)  =  cos  O  =  cos  (am(A,B)) 


sn(A,B)  =  sin  O  =  sin  (am(A,B)) 

At  this  point  we  get  an  expression  showing  the  character 
of  the  wave. 


(3-95) 


z  = 


H 


(1-  ek2)  +  (^k^)cn^('a]||,] 
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We  now  Introduce  two  i?estrlctlons  on  the  constants  Involved 
as  we  did  In  Case  1.   A  complete  cycle  of  (Z)  is  characterized  by 
a  cycle  for  (D)  equal  to  {v   radians).   The  corresponding  value 
for  the  cycle  of  (a)  will  be  called  (a). 


(3-96) 


w 


dO 


0  y  1  -  k^  sln^  O 


-'f^^ 


The  values  of  a  are  restricted  because  of  continuity  to  the 
following 


(3-97) 


77   2Tr 


(j)  =  1,2,5, 


where  (J)  is  the  number  of  crests  of  the  wave  in  the  trough. 


(3-98) 
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IT   ihi:  _ 


K 


We  combine  these  results  with  equations  (3-75,  95) 


(3-99)   e^i  =  eE^ 


2 
1  _  7  c 


p 


2HS^=n^(a/H|,; 


Setting  the  volume  of  fluid  with  the  wave  present  equal  to  the 
steady  state  fluid  volume  Is  equivalent  to  (3-100). 


(3-100) 


//    ET)-    d  Vol  =   0  =    //    eri-pdpda 


The   next   two   relations  will   be   required. 


(3-101) 


a 

/  ^"  =  "  -  *  i^ 


K 
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a 


(5-102)         /   cn^(Aa,k)da  =  j 


0 


•\< 


K 


K  ^  E 
k    k 


Equations  (3-99,  100,  101,  102)  now  give  us  (3-103) 


(3-103) 


eE, 


7H 


b^  +  2ab      2-^ 
-  7  c 


21n(l+b/a) 


K  - (l-^)K  -  ^E  =  0 


In  general  we  will  specify  the  height  to  depth  ratio  (qo/Qn ) 
and  the  number  of  crests  (J)  for  the  wave.   This  will  determine 
(erj^)  uniquely  except  for  a  special  case  discussed  immediately 
after  the  general  case.   Other  pairs  beside  (qo/q-]>(j))  could  be 
chosen,  but  this  will  not  be  discussed  here.   Note  the  limits  of 
the  range  of  (qg/q-^)* 


(3-104) 


0  <    (qg/q^)  ^  ^ 


With  (qp/q-i  )  specified,  equation  (3-85)  implies  that  (k  ,^) 
are  known.   With  (j)  specified,  equation  (3-9^)  .  piles  (H),  the 
amplitude,  is  known.   Now  equation  (3-103)  can  b^  used  to  determine 
(eEp)  except  for  this  special  case.   The  special  case  is 
characterized  by  (3-107),  a  special  value  for  c. 


(3-105) 


c 


c^  = 


b^  +  2ab 


27   ln(l+b/a) 


Equation  (3-103)  reduces  to  (3-106)  or  (3-10?) 


(3-106) 
(3-107) 


r\j  i\. 


(1-  ^)K+  ^E  =  0 


2     '\j        <v 

(k  -  2)K  +  3E  =  0 
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For  this  special  case  we  solve  (5-10?)  for  (k  )  and  hence 
for  (qp/Q-i  )  and  (4)»   In  this  case  we  specify  the  number  of  crests 
(J)  and  the  amplitude  (eE^);  then  (et)^)  Is  uniquely  determined. 

We  are  now  able  to  tabulate  some  of  the  dimensioned  results 
determined  up  to  this  point. 


(3-108)  V^  = 
(3-109)  V^  = 
(3-110)     V^  = 


_j|(H^)^/^  \  sin  20  ^1  -k^   sin^  O  i'    +0(e5/^) 

o  

_^(H^)^/2  i^  sin  20  Jl  -k^   sin^  O  ^g  +  0(e5/2) 

^   (7 +H(1  -  ^k^)  +  eE^/S +H^k^   cos^O 


.    (Hgk)^y 

■^  — znr- 


+   £^2   +  °(^^) 


5 -   2cos   aD(l  -  k      sin  O) 


(3-111)  T)      =     T^Q     +     ETl 


0(£^) 


4^eEil-4£)-ZHi-iKi) 
2p  \         p      /  p 


.  H2|!£  en^   (^,  R)   +.0(e2)    . 


(3-112)        P  =  6g(nQ  +  ET)^  +  e^ng-s)   +  O(e') 

I-    2p  \      p    y  p 


^^cn^    (^,k) 


+  O(e^)    . 
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Chapter  IV 

IRROTATIONAL  PROGRESSING  CNOIDAL  WAVES 
(CONSTANT  A/ELOCITY) 

In  Chapter  III  we  set  the  time  derivatives  equal  to  zero 
and  found  that  at  a  critical  condition  for  the  velocity,  waves  of 
cnoidal  shape  were  solutions  to  equations  (2-28,  33,  ^h ,    35,  36). 
These  waves  were,  of  course,  stationary  in  space.   In  this  chapter 
we  will  look  at  the  case  of  waves  rotating  in  the  same  circular 
trough  with  constant  velocity.   This  motion  is  characterized  by 
the  dependent  variables  being  functions  of  (p,s,a-Ot)  in  the 
original  coordinate  system,  or  (p,s,0-ooT)  in  the  stretched  co- 
ordinate system.   We  note  that 

In  terms  of  the  derivatives  of  the  dependent  functions,  this  wave 
may  be  characterized  by  (4-2). 

(k-2)  ^J=_S^J 

yf  n  =  -  »3Fi 

Under  this  transformation,  equations  (2-28,  33,  34,  35,  36) 
go  into  the  following  equations,  and  because  of  (4-2),  we  get 
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(E  =  constant).   The  problem  Is  to  solve  (4-3,  k,    5,  6)  and  then 
ro*^    ♦"he  pressure  by  substituting  the  solution  into  (4-7). 


(4-3) 

On  s  =  n(p.^)» 
(4-4) 

On  s  =  Ti(p,e), 
(4-5)        E 

On  s  =  G(p), 
(4-6) 


b^     ^  1    ^     ^    h    ^    e     b^' 


$  =  0 


%   -5p 


n„  =  ETl, 


-^  -  00 
.P 


=  e 


n.$,)|5 


+ 


![%''- (Is ''■ 


Is  -  lp°'  =  0  • 


(4-7)      TT  =  |[e  -|($p)^  -|(53)2]-s-l^.' 


$ 


X0 


2p 


^  -  OD 


We  now  make  the  same  substitutions  as  before  in  (3-11,  12,  I3) 
with  this  additional  note.   We  will  solve  for  the  dependent  variables 
as  functions  of  (p,s,9)  with  the  understanding  that  (Q-oot)  will 
be  substituted  for  (0)  when  we  get  the  solutions  at  the  end  of  this 
chapter. 

It  is  helpful  at  this  point  to  look  ahead  and  define  a 
particular  cross  section  of  the  fluid  in  the  circular  channel. 
Diagram  4-1  together  with  equation  (4-30)  establish  the  region  (  (f^.) 
of  undisturbed  fluid  at  the  critical  velocity.   Given  the  curve 
(s-G  =  0)  and  the  particular  value  of  velocity  (cu)  we  will  have 
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to  solve  for  (a,b,7).   The  upper  curve  will  be  seen  to  be  the 


leading  term  in  the  expansion  for  [t])  ;    it  will  be  named  (-q  )  later 
in  this  chapter. 

Let  us  first  look  at  the  (e  )-equatlons. 


(4-8) 


'4p^  +  ^p^)5d  =  °  -« 


(^-9)     1^   =0  on  h6l    [s  -T]      =  0,   and   s-  G  =  0} 

N  ^ 


(^-10) 


2Eo  =  4J'  +  (L)' 


Using  the  maximum  modulus  theorem  for  elliptic  equations  we 


get 


(^-11) 


$n  =  Fn(S)  . 


^  "  '0 


(^-12) 


Eq  =  0  . 


We  proceed  to  the  (e  )  equations. 


(^-13) 


^    0     ^     +    ^     D     ^' 


0 


ll-    -^     on£ 


On  s  =  G(p), 


(4-14) 

On  s  =  tIq, 


5i  ,.= 


0  . 


N 


(4-15) 
and 


5i  -  Ii  ^0 
s    p   p 


n 


0. 


F 
0 

~? 

.P 


03 
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(4-16) 


^0  =  %  -^  ^0 


CO  -  — 


0 


2p  J 


the  last  two  equations  we  get  on  s  =  t^  , 

■^2 


(4-17) 


s    p  p 


^0 


Applying  the  divergence  theorem  to  equations  (4-13,  14,  1?) 
we  find  that  'i     =   0.   Let 


(4-18) 


Fq  =  70  . 


This  implies  the  following  results 


(4-19) 

and 

(4-20) 


$1  =  F^(e)  , 


Tlo  =  E^  +  7 


Qi    - 


7 


2p"J 


Proceeding  to  the  (e  ) -equations  we  get 


(4-21) 
On  s  =  G, 


^^  P  ^  -^  ^  P  ^% 


(4-22) 


On  s  =  t)q. 


(4-25) 


and 


I. 


0 


■N 


Ip  -$5  1 

P  "P 


2.''0_  "  "^lo^T^ 


•0 


co)  , 
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(4-24) 


7 


P 


The  last  two  equations  yield  the  boundary  condition  (4-25). 
On  s  =  tiq. 


(4-25) 


$2  -5; 


s    p  p  p 


We  note  that  the  solution  of  equations  (4-21,  22,  25)  can 
be  written  as  follows: 


(4-26) 


f^{p,e,s)  =  F\{e)^{p,s)+F^i9)   , 


where  (Fp)  is  as  yet  undetermined  and  (^(p,s))  is  a  solution  of 
the  following  equations.   On  iR,  , 


(4-27) 


S    S  ^  S  ^ 


Y  =  -  - 
P 


On  s  =  G, 


(4-28) 


\  =  °  • 


On  s  =  TJQ, 


(4-29) 


s    p  '0 


=  -  (CD 


2^)' 


P 


To  determine  (F, )  exactly,  we  will  have  to  solve  this  last 
set  of  equations  for  the  particular  channel  (s-  G  =  0).   However 
we  can  still  determine  the  nature  of  (F-,  )  without  exact  knowledge 
of  the  solution  1. 
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divergence  theorem  to  (4-21,  22,  25)  we  get 
.....liar  blf..:  -   "^n  equation  (4-50);  this  equatl-..  ,^ 
•ous  to  (3-34). 

a+b 


(4-30) 


ff'-^-f         P(^-^)' 


dp 


<e 


=  0 


We  set  F,  ^  0  to  get  the  non-trivial  solutions.   This  gives 
us  an  equation  to  determine  the  values  of  (7)  for  each  channel 
(s-G  =  0)  and  for  each  angular  velocity  (00).   Integrating  with 
respect  to  (s)  in  (4-30)  results  in  (4-31)  which  we  then  integrate 
with  respect  to  (p)  resulting  in  a  quadratic  equation  (4-32  or  33 
or  34)  for  the  critical  values  of  (7). 


(4-31) 


/ 


a+b 


37CO  .   ~2 


dp  =  0 


(4-32) 


(f)^-  (y) 


31n(l+b/a)' 


ab(l+b/2a)_ 


2a'^(l+b/a)^ 


a+b 


G-E, 


+ 


7  ab(l+b/2a) 


dp  =  0 


TU  ^^^   f^f     /7\f2a^(l+b/a)^ln(l+b/a).    2a^(l+b/a)' 
^  ^^^^   (si  ~^Z)   b(l+b/2a) ^  ^  3 


_^  2a^(l+b/a)^ 
3co^b(l+b/2a) 


a+b 


G-E, 


dp  =  0 


a 
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Ik   ^lil      iy\  _a^(H-b/a)^ln(H-b/a)   +    |a°(n-b/a)''ln^(l+b/a)     2a''(l+b/a)^ 
'*-'*)      '^'=  bd+b/Sa)  -I  b^,i,b/2a)2 S 


CD 


3cn^b{l+b/2a)     ^  P 

3. 


With  the  same  understanding  as  in  Chapter  III  we  can  set 
(E,  =  0).   Let  us  briefly  examine  some  examples  of  the  solutions 
of  (4-32,  33,  34). 


Case  1: 


If  OD 


=  0,  we  get 


(4-35) 


7 


_  + 


'   ^       Q   a+b 
2a^(l+b/a)^   r    -G 
3b(l+b/2a) 


dp 


This  is  the  same  as  the  value  in  Chapter  III  equation  (3-38). 

Case  2; 

If  we  set  (7)  equal  to  the  value  in  (4-35)  we  get  both  the 
stationary  wave  (co  =  0)  and  a  progressing  wave 


(4-36) 


~  ^  3yln(l+b/a) 
ab(l+b/2a)   ' 


Case  3: 

Another  ca.se   of   interest   is    (7  =   0).      In   this   case   we  get 
equation    (4-37). 


(4-37) 


CD   =   + 


a+b 


-f^p 


a 


4   ab(l+b/2a) 
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It  is  also  interesting  to  look  at  the  solutions  as  (a)  gets 
large  with  (b)  relatively  fixed.  For  this  case,  let  (Ht)  be  the 
average  height  of  the  fluid.  The  result  is  relation  (4-58). 


(4-33) 


(au)  -  r/a)  —  (Ht)  as   (a)T  . 


The  term  on  the  left  is  the  square  of  the  relative  velocity  of  the 
fluid  motion  (steady  state)  to  the  wave  velocity  (  o  lowest  non- 
trivial  order  of  (e)).   The  relationship  {k-J>S)    in  the  limit  is 
the  same  as  the  critical  velocity  expression  for  infinite  straight 
channels. 

Returning  to  the  original  problem,  we  will  determine  the 
nature  of  (F-,(0))  and,  as  a  result,  we  will  determine  (Ti,(p,9)) 
from  the  (e  )-equatlons  which  follow.   On  (R. , 


On  s  =  G, 

(4-40) 

On  s  =  t]q, 

(4-41) 

On  s  =  tIq, 

(4-42) 


t 


=  0 


^N 


=  E- 


Ft^/2p^  +  $p  (S 


^-) 


P 


h  -  h  ^0 

L   3      P   PJ 


+  T- 


$2  -h     ^0 
ss   ^ps  p_ 


-h 


2F^r 


-   (03  -  ^)    I 


2F,F,   ^ 
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This  last  equation  Is  more  useful  after  the  substitution  of  (4-24, 
26) .   On  s  =  i]    , 


(4-43) 


i3  -^3^0   ='^'i 


-^(co 


2^) 


p 


F2(cD 


2      r*  ^  7' 
p       C  p^ 


1 


ss 


+  F  F 
■^11 


(S  -  -^) 


P   \P 


7^Y   \   27^  ■ 

¥     +       P^    +  __a 

ss     3   /    ^3 


The  divergence  theorem  Implies  (4-44). 


(4-44)  F^ 


a+b 


YdA 


^(OD  -  ^) 


R 


pdp 


+  F, 


a+b 

dA    r    /~    7 
-  /     (en  -  ' 

P    J 


2 


(to ^ 


)   pdp 


I-  je  a 

a+b  , 


+  F^F^ 


y   "¥      \       27^ 

~^M  ^   ss    3   /    3 
P   \P        P   /    P 


(^--4)  f4.-^__  +  — ^  )  + 


pdp 


s=ri 


0 


a+b  .^2^ 
+  F.E^  /    (  g^  -  pY 


1^2 


ss 


a 


dp  =  0 


s=ii 


The  coefficient  of  (Fg)  In  (4-44)  Is  zero  because  of  (4-30).   The 
coefficient  of  F**  is  simplified  by  (4-45). 
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JJ       P 


-r 


a+b 


P 


S=T1, 


pdp  =  /|  -^(V.pV4')dA  +'Pp'i'i^^ 


-ffi-n^' 


pW)  +V.  (p^W))dA 


Let 


r '  \ 


and 


=  r|^p{VT.W)dA 

=  /f  p((^3)^+(^p)^}dA  >  0 


(4-U6)    niQ  =Jfp[{^^)^   +    (^  )2]dA  , 


a+b 


(4_U7)    i\  =/   [(^  -   ^M  P^ 


2   ^'^ps 


ss   p 


m^ 


=  E, 


a+b  . 

a      ^ 


^2vj, 


ss    p2 


£_S 


dp  , 


S=T1, 


2y^ 


dp  , 


S=T1 


0 


(4-49) 


Y(e)  =  F^(e)  . 


With  these  definitions  for  {mQ,m,  jin^jY) ,  equation  (4-44) 
reduces  to  (4-50) which  is  the  basic  equation  for  cnoidal  waves 

(3-65). 


(4-50)  or  O63) 


itIqY  =  m^YY  +  mgY 


The  discussion  of  cnoidal  waves  following  equation  (3-63)  applies 
here. 

For  the  special  case  7  =  0  we  have  the  following: 
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(4- 

-51) 

(^■ 

-52) 

(4- 

-55) 

(4- 

-54) 

i^- 

-55) 

(4-56) 


Tlo   =   0 


•n^   =  Eg   -f  Yo) 


ri„   =  E^   -  ^— ^  +  toY^   +  oiF, 
2  5        2p^ 


m. 


0 


=irp[<*si^  +  (*p)^]dA 


m-,    = 


/e 


03 


a+b 


[pY        -  -] 
•■^    ss        p -■ 


dp 


a-fb 
a 


s=no 
dp 


S=T1 


0 


(4-57) 


itIqY  =  m-|_YY  +  m^Y 


Recapitulating   the   results   for   this   chapter,    we  have 


(4-58) 


'§_=ye+el    Yde  +  e^[^Y  +  F2]  +0(e^)    , 


(4-59) 


T]     = 


E,    -t-  7CD   -  ^ 


^P^J 


+   e 


Ep  +  Y(33  -  ^) 


+   e 


-  ^  +    (o)   -  ^)(YY+F^) 


^        2p^ 


P 


+   O(e^)    . 
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Chapter  V 
LIMITING  CASE:   STRAIGHT  CHANNEL 

At  the  beginning  of  this  report  we  noted  the  lack  of  an 
existence  theorem.   In  order  to  establish  some  degree  of  credence 
to  this  method,  we  shall  show  that  In  the  limit  these  solutions 
go  smoothly  into  the  solutions  for  straight  channels.   Since  we 
have  an  existence  theorem  [5]  for  the  solution  in  the  straight 
channel,  we  might  expect  existence  for  the  waves  discussed  in  this 
report.   We  will  show  that  as  a  -*  00 ,  b  fixed,  the  solution  to  the 
irrotational  motion  goes  into  the  solution  discussed  by  A.  S. 
Peters  [2]. 

The  general  solution  in  [2]  which  is  described  with  respect 
to  an  Inertlal  system  traveling  with  the  wave  is  specified  as 
having  the  flow  at  infinity  equal  to  zero.   This  is  equivalent  to 
the  case  in  Chapter  4  where  7=0.   This  Is  described  by  equations 
(4-37,  51  to  57). 


(4-57) 


2 

00 


-     a+b  1    /  r     ^"^^ 

La  J  /    l-  a 


pdp 


(4-51) 


^n  -. 


0  =  0 


(4-52) 
(4-53) 


Ti^   =   Eg  +  cuY 


0  O  '^^•  ~' 

Ti^  =  E,   -   YV2p      +  cdY^   +  cuF, 
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(4-54)  m^  = 


(4-55) 


0 


viu    = 


a+b 


CD 


f  (P^  -    -) 

J  ^^    ss        p' 


dp   =  a3[pY   (p,G)  +  21np 

r 


'0 


a+b 


a+b 


(4-56)  m^  =  E_    /  p^ 


2  2 


ss 


dp    =  E[p^    (p,0) +  Inp 

r 


Tlr 


a+b 


(4-57) 


m„Y   =   m-,  YY   +  iiigY 


The  auxiliary  function  (Y)  is  defined  by 


(4-27) 
(4-28) 

(4-29) 


1  a      a  ^  s^ 


Y  =  - 


1 
p 


n  oe 


Yj^  =  0     on   s  =  G 


^N  =  """    on   s  =  kiq 


First  let  us  look  at  the  wave  speed.   The  linear  velocity 
corresponding  to  (oo)  is  (aco).   In  the  limit  as  a  f  oo  we  get 


a+b 


(5-1)   (acD)' 


F 


-G  =  the  average  height  of  undisturbed  fluid. 


In  [2],  A.  S.  Peters  uses  a  second  parameter  to  reduce  the  average 
height  to  unity.   This  corresponds  to 


(5-2) 


(acjo)  -»•  1  . 


Let  us  now  substitute  (4-52)  into  (4-54)  to  (4-57).   This 
gives  us 
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(5-5) 
(5-4) 
(5-5) 
(5-6) 


^0^1  =  n^]^^!^!  +  f"2^i  . 


Q=^{/^p[i^2^y  2]dA  , 


m. 


^  =  (p*  (p,0)  +  31np) 


a+b 


m^  =  2E2ln(l +  b/a)  . 


Expressing  (9)    in  terms  of  (a9),  the  linear  dimension,  in 
(5-3)  gives  us  (5-7,  8) 


(5-7) 


(5-8) 


^3 
/K    a 
m, 


ii    d 


0  ^3^^i  =  "^^1  dlley  ^1  ^  "^2  duey  ^i 


m_^  =  a  m^  ,   ^i  ~  ^^  '    ^p  ~  ^^p  • 


Let  me  now  substitute  T*  =  -a  "¥   in  (4-27,  28,  29). 


(5-9) 


1  s  ^  a' 


L^P 


?  +  p  ^  + 


Ss' 


m*       a 


in 


(5-10) 
(5-11) 


n 


r    ~\2 

(aa>)    on  r]      , 


^*  =  0 
n 


on   s  =  G  . 


Taking  the  limit  as  a  -►  co  we  get  equations  identical  to 
those  on  [2],  page  462  for  (^*). 

/V       <«>       /N 

Now  let  us  express  m^,    f?L,  A-  in  terms  of  {"H*)   and  then  let 

a  — »  J-  . 
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(5-12)    \  =  ff  {p/s^)[^f  +  ^f]dpds^ff[^f  +  ^f]apdi 


(I  <l- 

a+b 


-    3aln(l+b/a)  -*  -3b   +   [Tp(a+b)    -  Y*(a)] 


(5-13)      m^   =    (p/a)^^ 

(5-1^)      trig  =  2aE2ln(l   +  b/a)  -y  abE^ 

(m. )  corresponds  to  (m. )  in  [2],  pages  456,  457,  the  Irrotational 
case.   We  now  see  that  (5-7)  is  identical  to  [2],  eq.  3.47.   This 
equation  expresses  the  total  family  of  cnoidal  and  solitary  waves. 
The  limiting  processes  used  are  valid  since  we  are  dealing 
successively  with  uniformly  bounded  continuous  functions,  elliptic 
P.D.E.  on  a  smoothly  changing  family  of  domains,  and  an  O.D.E. 
with  coefficients  which  are  uniformly  bounded  above  and  below  for 

the  unbounded  parameter  (a). 

(a) 

Since  the  O.D.E.  for  r),^   (a9)  goes  smoothly  into  the  equation 

for  Ti|   ^(x),  a  simple  examination  of  the  solutions  in  terms  of 
the  coefficients  shows 

(a)  /  a\  (cD  )  /   ^ 

r|^  '{B.Q)  -*  T]^      '(x)  . 

Higher  order  terms  may  be  handled  in  a  similar  manner,  and  we  may 
also  treat  the  rotational  case  this  way. 
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Chapter  VI 
ROTATIONAL  CNOIDAL  WAVES 

This  class  of  problems  Is  handled  the  same  as  In  Chapter  IV. 
We  assume  that  the  dependent  variables  are  functions  of  (p,s,a-nt) 
or  (p,s,9-a)T)  in  the  original  or  stretched  coordinates  respec- 
tively. 

We  note  again  that 


(6-1) 


<^  =  0//g  . 


Equations  (2-19)  thru  (2-25)  imply  (6-2)  thru  (6-8).   In  the 
region  of  fluid  we  have 


(6-2) 


^(pU)  +^(PW)  =  -e   ^(V)  , 


^ 


(6-3) 

(6-4) 

and 

(6-5) 


UU   +WU   +eTT   -  —  +  Uo(-  -  o)) 

p    s    L  p   P    ^  p 

U  ^(pV)  +  WV  +  e  — ^  +  Vfl(-  -  Z) 
p  dp  ^        S    L  p       9  p 


UW   +  WW   +  e 

P       3 


[i-K 


=   0  , 


=  0  , 


=  0  . 


On  (s-T)(p,9)  =  0)  we  have 


(6-6) 

and 

(6-7) 


Tr=  0 


W  -Ut]^  = 


-  (d) 
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On  (S-G(p)  =  0)  we  have 

(6-8)  W-  UG'  (p)  =  0  . 

If  G  (p)  =  00,  i.e.,  vertical  walls,  we  must  set  (U  =  0)  at  these 
walls  in  addition  to  (6-8)  elsewhere. 

We  make  the  substitutions  given  by  (6-9)  and  proceed  to  solve 
the  resulting  equations  recursively  in  increasing  order  (i)  of  (e"""), 

00 

(6-9) 


u  =  y^  U.  e^ 

00  . 


1=0 

00 


i=0      ^ 


00 


i^   ^ 

CO  . 

TT-ZZTTiB^ 

i=0        ^ 


If  we  try  to  proceed  in  all  possible  generality,  we  find  that 
the  solutions  fall  into  many  categories  that  must  be  considered 
separately.   The  steady  state  (low  order)  condition  that  is  being 
investigated  causes  this  breakdown. 

In  all  the  cases  that  we  have  considered  and  I  suppose  in  all 
possible  cases  there  are  three  familiar  results.   The  first  is  the 
bifurcation  in  the  equations.   This  occurs  consistently  showing  a 
dependence  on  a  critical  speed.   The  last  result  is  the  cnoidal 
shape  of  the  waves  in  the  direction  of  the  angular  variable.   In 
most  cases  there  is  a  horrendous  bookkeeping  problem  with  lengthy 
equations. 
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Instead  of  looking  at  all  these  results,  let  us  Just  look 
at  one  Interesting  simple  case.   The  simple  nature  of  this  case 
is  by  no  means  representative  of  this  class  of  solutions,  but  it 
does  Illustrate  the  surface  waves  that  we  are  interested  in. 

Let  us  assume  a  rectangular  channel  (see  Dlag.  6-1), 
stationary  waves  (a>  =  0)  and  a  particular  steady  state. 

Assume: 

(6-10)  G  =  -h  <  0  ,   for  a  <^  p  ^  a+b  , 

riQ(a)  =  0  , 
CO  =  0  , 
Xq   =  /hp/a  , 
Uq  =  Wq  =  U^  =  W^  =  0  . 

Let  us  note  that  the  assumptions  in  (6-10)  are  not  contradictory 
in  any  way,  but  describe  a  particular  situation. 
The  ( e  )  and  ( e  )  equations  imply 

(6-11)  TTo  =  s  =  tiq  =  ^(p-a) 


which  is  an  exact  steady  state  solution.   We  also  note  that 
^0 


2 
(^1^  -  G  =  Xn  )  ^°^  this  case;  the  assumptions  were  set  up  so  that 


p 
this  critical  condition  would  be  identically  satisfied.   The  (£  ) 

equations  are  now 


(6-12)       ^(pU^)  +  ^(PW^)  =  -Xi 


d 
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(6-13) 

^P^ 

2X0X1 
p 

:  ax^yh/pa 

{6-n) 

^s   ^ 

0  , 

(6-15) 

'ie  = 

->^O>^l0 

3U. 

(6-16) 

TTi  -  Til 

(6-17) 

Wg-UgTlc 

(6-18) 

Wg  =  0  , 

(6-19) 

Up   =   0    , 

=  -Xt  yhp/a 5- j/hp/a    ,        in  region   >.^  , 


on    s   =  tIq    , 


UpTlQ      =  Tl^   Xq/p   =  '^i  /h/pa    ,    on   s   =  i1q   , 
p  0  0 


on   s   =  G  =    -h   , 


on  p   =  a,    p   =  a+b 


Equations  (6-12)  thru  (6-I9)  have  the  following  solution  (6-20) 
thru  (6-23),  in  terms  of  functions  (A(0))  and  R^(p).   The  first 
of  these  will  be  determined  from  higher  order  considerations.   The 
second  is  arbitrary  and  will  not  effect  the  class  of  solutions  for 
(A(9)),  but  will  determine  the  specific  solution. 

(6-20)  TTi  =  Ti^(p,0)   =  A(|l  +R    (p) 

p 

(6-21)  Xi-p^Ti      =  I  Ti    /p^7h  =  |/pi7h  [r!  -  % 

^0        p  p  ^        p-^ 

(6-22)  Up   =   0 

(6-23)  Wp  =    -|^    A{e)  = —^ — A(0) 
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The  solvability  of  the  above  problem  (restated  in  (6-24)) 


2 
was  assured  by  the  fact  that  (t^,.  -  G  =  Xn  ) »  ^^®  critical  condition 


for  this  case. 

u   w 

(6-24)  V  (T^,r^)  =  f(p,0,s)   in  R.   , 

plus  boundary  conditions. 

The  {^    )   equations  now  show  us  the  cnoldal  character  of  A^d). 

(6-25)   ^(pU^)  +  ^(pW^)  =  -X2^  , 

(6-26)  TT^  =  -(^)X  , 

s      p 

(6-27)  TT^    =  ^l!l^  =  2x2/h7p^  +  1^  fR'  -  ^1'  ' 

P  ^  9 

(6-28)  TT^^  =  -XoX2^  -  XiXi^  -  u^^(pxo) 

=  -X2^/Hp7i:  -  I  Uy^h^  -   fn  ^1  ^1  g 

=  -X2  /HpT^  -  I  U/Hp7a  +  -^  A  [r'  -^^\    ,      in  ^ 

(6-29)   172  =  ^2'  on  3  =  tIq 

(6-30)   W^  -Ti^^U^  =  -^1^2/^20  Vp  +  ^l0  Vp 

2AA  /aTh  _^  /.  /^  ^ 
P  ^ 

^iy^  tR72  _  R/p]    on  s  =  ^0 
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The  method  we  use  now  Is  to  integrate  (6-26)  to  solve  for 
(TTp)'   W^  then  use  the  divergence  theorem  to  get  a  differential 
equation  for  A(0). 


(6-31)    J\o   =  no(p.^)  -  \ 


2     2 


+  h(s  -tIq) 


=  Tip(p,e)  -  ^  [s2+2hs  +  h^(l-  (p/a)2)]  . 


We  see  that 


(6-52) 


-"^O  '  ''^O 


=  V 


pu^  pW 


3 


PXo  '  PXo^ 


U- 


^0 
X2 


PXo 


PXn   77^ 


'0    PX 


TTq   +  X2  v/hp/£ 


-e        ^e 


+  — ^  A  L-^  -  R 

2hp      p 


■^] 


We  now  apply  the  divergence  theorem  to  (6-33^  3^)  to  get  (6-35) 


'U,  W,\       ^2     „. 
(6-33)    V.(-^   5=  ^-^a){^.   R') 
''^•0   -^0  /     2h^p^     p-^ 


+  -|-U    -  ^(s^+2hs  +  h^(l-  p/a)^)i 
p  h  (.  ^0   2p^  J 
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(6-34)    W,  =     on  s  =  -h 


U,  =  0   on  p  =  a,  p  =  a+b  . 


W 


!S   _   "3  _   2AA(a/h)    1       A(a/h),  ■,     /. 


'p  ^0 


p  '-=^0 


(6-35) 


on  s  =  T)-. 


dp 


S=T) 


0 


Evaluating  the  first  integral  we  get 


'^-'^'  /^•(S'S) 


r  -^  A  {^  -  R')dp+  Tcti  /p)dp 

^a     2hp^         p^  ^^       ^0 

b       ...        > 
^     ^\    /    [s^+ 2hs+h^(l  -  p/a)^]dsdp 


2hp 


-h 


Integrating  with  respect  to  s'and  p  in  equations  (6-35,  36)  gives 
us  (6-37,  38). 


a+b 


(6-57)  O.a^A-/   S^ 


^I8a^r    d^_6a^/^   J^^(R7p)dp 


P'     h-   "^„    p-    h-  -     P 


3  J      ;?v 

a 


(6-38)   A 


1  - 


(l+b/a)^J 


18  AA 


(l+b/a)^_ 


+ 


^,r7>'/p)ap 


Equation  (6-38)  is  the  form  of  (3-63)  with  m,  negative,  hence 
the  cnoidal  nature  of  A(9).   This  is  precisely  the  result  that  we 
want. 
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Diagram  2-1 
Coordinate  System 


T. 


a,9,T^'    T^ 


P(p,a,s) 
P(p,9,s) 


a+b 


fixed  surface:   s-G(p)=0 
fluid 
free  surface:   s-Ti(p,9)=0 


Diagram  3-1 
Definition  of  Region  ^ 


fixed  surface:   s-G(p)=0 
fluid  (steady  state)  region 


free  surface 


Note: 


The  value  of  r]      =   ^^^  is  in  anticipation  of  equation  5.25. 
It  is  presented  here  so  as  not  to  obstruct  the  procedures 
later. 
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Diagram  3-2 
P(Z) 


p(z)=  xz' 


(q2-<^i) 


(-Qo) 


Region 

of 

Interest 

P(Z)>0 


(Ql) 
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Diagram  4-1 

Definition  of  Region  (l 


a 


a+b 


■*^  P 


o 

•H 
-P 

:3 

rH 
O 
> 

u 

o 
m 

•H 


fixed  surface:   s-G(p)=0 


fluid  (steady  state)  region 


2 


free    surface:         s-ti„  =  s  -  E-,  -  cu7+ -^=  0 

P 


Note:   The  expression  for  (tIq)  Is  In  anticipation  of  equation  4-20. 
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Diagram  6-1 
Region  /2 
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